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ABSTRACT 

"Hard" massive black hole (MBH) binaries embedded in steep stellar cusps can shrink via three- 
body slingshot interactions. We show that this process will inevitably be accompanied by a burst of 
stellar tidal disruptions, at a rate that can be several orders of magnitude larger than that appropriate 
for a single MBH. Our numerical scattering experiments reveal that: 1) a significant fraction of stars 
initially bound to the primary hole are scattered into its tidal disruption loss cone by gravitational 
interactions with the secondary hole, an enhancement effect that is more pronounced for very unequal- 
mass binaries; 2) about 25% (40%) of all strongly interacting stars are tidally disrupted by a MBH 
binary of mass ratio q = 1/81 {q — 1/243) and eccentricity 0.1; and 3) two mechanisms dominate 
the fueling of the tidal disruption loss cone, a Kozai non-resonant interaction that causes the secular 
evolution of the stellar angular momentum in the field of the binary, and the effect of close encounters 
with the secondary hole that change the stellar orbital parameters in a chaotic way. For a hard MBH 
binary of 10^ Mq and mass ratio 10~^, embedded in an isothermal stellar cusp of velocity dispersion 
(T* = 100 kms~^, the tidal disruption rate can be as large as iV* ~ 1 yr"-'^. This is 4 orders of 
magnitude higher than estimated for a single MBH fed by two-body relaxation. When applied to 
the case of a putative intermediate-mass black hole inspiraling onto Sgr A*, our results predict tidal 
disruption rates ~ 0.05 — 0.1 yr~^. 

Subject headings: black hole physics - methods: numerical - stellar dynamics 



1. INTRODUCTION 

Close massive black hole (MBH) binaries are expected 
to form in large numbers followi ng the hierarchical as- 
sembly of massive ga laxies (e.g. iBegelman et al.l 119801 : 
IVolonteri et al. II2003D . but their merger history remains 
poorly understood. Few observational probes of the 
processes that lead to and accompany the shrinking 
and inspiral of a MBH binary have been proposed to 
date: 1) the gravitational slingshot ejection of hyperve- 
locity stars from the Galactic Cent e r into the halo (e.g . 
Yu & Trcmaine' '2003'; Twrn '2006': I Sesana et all [20061 . 
2007: Bro wn et a l. 2009: Perets 200"9| ): 2) interruption 
or redirection of jets due to MB H binary-accretion disk 
interaction or M]3H coalescence (iMerritt fc Ekerd l2002l : 
iLiu et all [200l iLiul [2003 : iLiu fc ChenI I2007D : 3) the 
coalescence of MBH pairs with masses in the range 
(10^ — 10^)/(1 -I- z) Mq giving origin to gravitational 
wave events that are one of the primary targets for 
the planned Laser I nterferometer Space Antenna (LISA ; 
e.g. Ea chnclt 1994: lHughed[2?M W^he fc Loebl[2003l : 
ISesana et al.l |2004L 120051 ): 4) the electromagnetic after- 
glow from a circumbinary accretion disk that would fol- 
low such coalescence (e.g.lMi losavli e vie & Phinncv 200j: 
Dotti et al.ll2006t iLippai et al...2008 : .Shields fc Bonni^ 



2008f) : and 5) the high-velocity recoil experienced by 



the plunging binary due to the asymmetric emission of 
gravitational waves (e.g. iBaker et all l2008l and refer- 
ences therein). A recoiling hole that retains the inner 
parts of its accretion disk may have fuel for a long- 
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lasting luminous phase alon g its trajectory, and shine 
as an off-center AGN fe.g. iMadau fc QuataertI l2004l : 
iBlecha fc Loebl[2008l: IVolonteri fc Madaull2008ri . 

In this Letter we return to the dynamical processes that 
determine the decay of MBH binaries in a stellar back- 
ground, prior to the gravitational wave regime, and put 
forward another possible observational signature of close 
binaries in the nuclei of galaxies. Using results from scat- 
tering experiments we show that gravitational slingshot 
interactions between an unequal-mass "hard" binary and 
a bound stellar cusp will inevitably be accompanied by 
a burst of stellar tidal disruptions, at a rate that can be 
several orders of magnitude larger than that appropriate 
to a single MBH fed by two-body relaxation. The dura- 
tion of the phase of enhanced tidal disruption is of the 
order of 10^ - 10^ yr. 



2. SCATTERING EXPERIMENTS 

Ana lytical techniques have been used bv llvanov et al.l 
(|2005| ) to study the enhanced stellar disruption rates in- 
duced by the secular non-resonant interaction with a non- 
evolving MBH binary. Here, we perform detailed numer- 
ical experiments of the close encounters between stars 
and the pair of MBHs, collisions that perturb stellar or- 
bits in a chaotic way and scatter stars initially bound 
to the primary MBH into its tidal disruption loss cone. 
Consider a MBH binary of mass M = Mi + M2 = 
Mi{l + q) (A/2 ^ Ml) and semi-major axis a, orbit- 
ing in a background of stars of mass m,, radius i?*, 
and velocity dispersion cr,. When a < Oh = GM2/4crJ, 
the "hard" binary lo ses orb ital energy by three -body 
slings hot interactions (jOuinla n 1996; Sesana et al.! l2006l . 
l2007t ). For unequal-mass pairs, the radius of influence 
Hnf = G(Mi + A/2)/(2crJ) is much larger than the hard- 
ening radius a^, and almost all interacting (low angular- 
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momentum) stars are bound to Mi.^ In the case of an 
isothermal stellar distribution around Mi, the total stel- 
lar mass within ah is equal to M2/2. 

The integration of the three-body encounter equations 
is performed in a coordinate system centered at the lo- 
cation of Ml . Initially the binary (of mass ratio q) has 
eccentricity e and a randomly-oriented orbit with M2 at 
its pericenter. Stars initially move in the x — y plane with 
pericenters along the positive cc-axis and random orbital 
phases. The initial conditions of the restricted three- 
body problem problem are then completely defined by 6 
variables, 3 for the binary and 3 for the star: 1) the in- 
clination of the binary orbit, 9; 2) the longitude of M2's 
ascending node, /; 3) the argument of M2's pericenter, 
0; 4) the semi-major axis of the stellar orbit, a*; 5) the 
specific angular momentum of the star, j^,; and 6) the 
orbital phase of the star, p,. We start each scattering 
experiment by generating 6 random numbers, with cosO 
evenly sampled in the range [—1,1], and I and </> uni- 
formly distributed in the range [0, 27r]. We sample a» 
logarithmically around a in the range [l/2a, 2a] where 
three-body interactions are strongest, randomly be- 
tween and 1 (corresponding to an isotropic distribu- 
tion), and evenly between and 1. The equation 
of motion are integ rated using an exp licit Runge-Kutta 
method of order 8 (jHairer et al.1119871 ). with a fractional 
error per step in position and velocity set to 10~^^. We 
have tested our cod e by reproducing Figures 4 and 6 in 
iSesana et al.l (|2008( ) and found excellent agreement, and 
run 10^ scattering experiments for each binary configura- 
tion. During each experiment the minimum separation 
^"111111 between the star and the MBH pair is measured 
and stored: stars on orbits intersecting the tidal radius 
of hole M, (i = 1,2), 

will be tidally disrupted at pericenter passage (neglect- 
ing general relativistic effects that set in when Mi ^ 
10^ Mq, see lHiilsinJTSl) . The results of our numerical 
experiments are shown in Figure [T] for binaries with dif- 
ferent mass ratios and eccentricities, all at separation 
a = a/j. The fraction of interacting stars that are scat- 
tered by M2 to within a pericenter distance r,„i,i < rn 
from Ml and arc tidal disrupted can be, for very unequal 
mass binaries, orders of magnitude higher than the corre- 
sponding number were Mi not in a binary. The latter is 
simply given by all the bound stars within Mi "tidal loss 
cone" , the region in the (a, — j, ) phase-space bounded 

by 

.2_/l {a*<rti) 

^i-"\2(ni/a,)2(a,/rti-l/2) (a, > m), 

where j, is the specific angular momentum of the star 
normalized to the angular momentum of a circular or- 
bit with the same semi-major axis. For a binary with 
q = 1/81, e = 0.1, the probability that a close encounter 
with a star having l/2a < a* < 2a results in a tidal dis- 
ruption by Ml (r„ii,i < 3 x 10~^a « rti) is more than 2 

^ Note that, for extreme mass ratios g <C 1, the encounter is 
essentially a two-body scattering as the star and the secondary 
hole move in the static potential of the primary. 



orders of magnitude larger than if Mi were single. The 
figure also shows that: 1) many more stars are disrupted 
by Ml than by M2; 2) the disruption probability de- 
creases with increasing q. This is both because the ratio 
rti/ ah decreases with increasing g, and because, as the 
perturbing force of the secondary increases, more stars 
are ejected altogether rather than disrupted; and 3) more 
stars are scattered into Mi 's tidal loss cone with decreas- 
ing binary eccentricity. 

3. BASIC THEORY AND DISRUPTION RATES 

Figure [2] shows an example of a three-body scattering 
leading to a tidal disruption after many pericenter pas- 
sages. In the presence of a secondary black hole, strongly 
interacting (a* ~ a) stars in nearly circular, highly in- 
clined orbits relative to the binary orbital plane undergo 
a slow secular evolution that periodically increases their 
eccentricity (in exchange for a lower inclination) and 
eventually drives them w ithin the tidal disr uption loss 
cone of the primary hole (jivanov et al.ll2005l ). Our scat- 
tering experiments show that this mechanism - analo- 
gous to_the sallied "Kozai effect" of celestial mechan- 
ics (jKozai I [l962[ ) - contributes but does not dominate 
the fueling of the tidal disruption loss cone in the case 
of very unequal mass binaries. The stars supplied to the 
disruption loss cone by the Kozai effect arc all those hav- 
ing normal component of the angular momentum, jz*, 
within the loss cone, i.e. all those stars within a wedge- 
like region in phase space where \jz*\ < jic- For a given 
stellar semi-major axis, a, <^ a, the fraction of stars that 
lie outside the tidal disruption loss cone but inside the 
"Kozai wedge" is 

/if(rti,a*)= / [ = 2jic - 2jf^. (3) 

When rti ^ a*^ fx — 2\/2rti/a*, which is ^J2a^/rti 
times larger than the fraction of stars already in the tidal 
loss cone, jf^ ~ 2r(i/a,. This result explains why the 
probability of stellar disruption for bound stars is much 
higher in MBH binaries than in single black hole sys- 
tems. It also predicts that only a fraction 2^2rtil a^, = 
(0.05, 0.09, , 0.15) for q = (1/81, 2/243, 1/721) of all stars 
with a/2 < a^, < 2a will be supplied to the tidal loss 
cone by the Kozai effect. Figure [1] shows that for very 
unequal mass binaries the tidal disruption probability is 
much larger than the above estimate. This discrepancy 
highlights the importance of close, resonant encounters 
with the secondary hole, which change the stellar or- 
bital parameters in a chaotic way and fuel the tidal loss 
cone. Figure [2] (bottom right panel) depicts the initial 
distribution in the (a* — jz*) plane of all the stars that 
arc disrupted in our numerical experiments. It is clear 
that the majority of disrupted stars lie outside the Kozai 
wedge. 

We now show that the ejection of ambient stars by 
a hard MBH binary will be accompanied by a burst of 
stellar tidal disruptions, at a rate that may be orders 
of magnitude larger than that appropriate for a single 
MBH. In our numerical experiments the time when a 
star first crosses the tidal radius of the primary hole is 
stored and used to calculate a disruption frequency. To 
translate this number into a stellar disruption rate in 
physical units one needs to specify the parameters of the 
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Fig. 1. — Upper solid line: Closc-cncountcr probability for bound stars interacting with the primary member (of mass Mi) of a MBH 
binary of mass ratio q, eccentricity e, and separation a = ai^. The vertical axis shows the fraction of stars N/Ntot with closest approach 
distance r^i^ < r. Lower solid line: same for M2. Dotted line: same in the case of an isolated MBH of mass Mi. Stars are drawn from a 
spherical isotropic distribution bound to Mi and have semi-major axis in the range a/2 < a» < 2a. The short vertical linos mark the tidal 
radii of Mi and M2 and the Schwarzschild radius of Mi (all in units of a) for Afi/m« = lO'^, r» = Rq, and cr* = 100 kms^^. From the 
top left to the bottom right, the fraction of stars tidally disrupted (rmin < rti) by Mi is 0.24, 0.19, 0.39, and 0.48, respectively. 



MBH binary and its stellar cusp. At the hardening radius 
GM2 
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the orbital period of the binary is 



Ph = 2^r^ 



G{Mi + M2) 



(3.3x10^ yr) a^^^Mj 



1/2 



1 + 9 
(5) 

where ctioo = cr*/100 knis^^ and M^ = Mi/lO' M©. If 
we assume now, for simplicity, that the stars bound to Mi 
follow an isothermal distribution, p*(r) = al / {2TiGr'^) . 
the total stellar mass between a/j/2 and 2ah is = 
3gMi/4. Normalizing the interacting mass in our nu- 
merical experiment to the isothermal case and rescaling 
the stored disruption times according to equation ([5]) , we 
obtain the stellar disruption rates shown in Figure [3l 

Although the standard Kozai theory does not strictly 
apply to strongly interacting stars, we use it here to de- 
rive an analytical scaling for the stellar disruption rates. 
The Kozai timescale at au is approximately 

37rg \ahj 

(jlnnanen et all 119971 : [Kiseleva et al.1ll998f) . The stellar 
disruption rate can then be estimated as 



where fx is the fraction of stars in the Kozai wedge and 
A ~ 0.2 is a correction factor accounting for the uncer- 
tainty in Tk and for stars that are actually ejected before 
disruption, as well as for properly weighting our scatter- 
ing experiments for the case of an isothermal profile. The 
numbers provided by equation ^ arc in good agreement 
with the numerical rates. The tidal disruption plateau, 
however, lasts much longer than (indicated by the 
vertical ticks in Fig. [3]), because at late times the dis- 
ruption rate is dominated by chaotic scatterings. 

The tidal disruption rates wc compute are many orders 
of magnitude higher than those. 



(2 X 10"'' yr"^) 



-1/3 



R© 



1/4 



K 



(9) 

derived for a single M BH fed by two-body relaxation 
(jWang fc Merrittll200l) . Note that in our calculations 
we have only considered stars with a/2 < < 2a. Tak- 
ing into account stars in a larger range of semi-major 
axis would further boost the binary disruption rates. 

4. CONCLUSIONS 

We have used results from numerical scattering ex- 
periments and shown that the tidal disruption rate in 
a stellar cusp containing a 10'' M© MBH binary can be 
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Fig. 2. — Top and botton left panels: Example of a chaotic three-body scattering leading to a tidal disruption. A star on a bound orbit 
of semi-major axis a* = 1.2a and eccentricity e» = 0.5 interacts with a MBH binary of parameters Mi = lO'^ Mq, q = 1/81, e = 0.3, 
a = ai^. Top left: stellar (black) and M2 (green) trajectories in the x — y plane. The primary hole Mt is located at the origin. Top right: 
same projected onto the x — z plane. Bottom left: separation between the star and Mi as a function of time (in units of the binary period 
P). The dotted line marks the tidal radius rti- Bottom right panel: Bound stars in the (a, — jzt) plane that are tidally disrupted during 
the interaction with a MBH binary of parameters Mi = 10^ M©, q = 1/81, e = 0.1, a = af^. The solid lines mark the boundaries of the 
Kozai wedge \jzt \ < jic- 



as large as 1 yr"-*^ over a timescale of ^ 10^ yr. This 
is orders of magnitude larger than expected in the case 
of single MBHs. After a tidal disruption, about half 
of the debris will be spewed into eccentric bound or- 
bits and fall back onto the hole, giving rise to a bright 
UV/X-ray outburst that may last for a few years (e.g. 
lReeslll988D . "Tidal flares" from MBHs may have been 
observed in several nea rby inactive galaxies (jKomossa I 
120021: lEsquei et al.ll2007j ). The inferred stellar disruption 
frequency is ~ 10~^ yr ~^ per galaxy (with an order of 
magnitude uncertainty, iDonlev et al.l |2002[ ) . The much 
enhanced disruption rates we have found here for MBH 
binaries can then be used to constrain the abundance 
of clo se MBH pairs in nearby galaxy nuclei (jChen et al.l 
|2008[) . It is interesting to scale our results to the scatter- 
ing of stars bound to Sgr A* , the massive black hole in the 
Galactic Center, by a hypothetical insp iraling companion 
of int ermediate mass ()Yu fc Tremainei '2003: S esana et al.l 
|2007() . The stellar density profile around the Galactic 
Center can be described as a double po wer-law, with 
outer slope ~ — 2 and inner slope ~ — 1.5 (iSchodel et al.l 
|2007() . If the density profile inside the influence radius of 
Ml is shallower than isothermal, 

P*{r < Tinf) = p*{ri„[) [ — ] (10) 

Vine J 

with 7 < 2, then the stellar mass between 2afi and 0^/2 
decreases by a factor (a/j/rinf)^"'', and the stellar disrup- 



tion rate in equation ([5]) decreases by a factor {q/^Y 
relative to the isothermal case. Using Mi = 4 x 10^ M0, 
cr* = 100 kms-i, 7 = 1.5, and 1/243 <q< 1/81, yields 
rates in the range iV* ~ 0.05 — 0.1 yr~^. 

There are a number of uncertainties in our calculations 
that require clarification before a firm statement can be 
made on the rates and duration of stellar tidal disrup- 
tions expected in galaxy nuclei hosting MBH binaries, 
and on the constraints imposed by the very low level 
of activity observed in the Galactic Center. First and 
foremost, our estimates of the tidal disruption rate as- 
sume a fixed binary separation a. In reality, both binary 
separations and eccentricities will e volve due to three- 
body slingshots (jSesana et al.|[2008[ ). This changes the 
Kozai timescale of the system and replenishes the supp- 
ply of strongly inter acting stars. According to Figure 7 of 
iSesana et al.l (|2008[ ). the evolutionary timescale of a bi- 
nary with initial eccentricity 0.1 embedded in an isother- 
mal cusp is th ^ q'^/^P. At a = ah we derive 

th ~ (3.3 X 10^ yr) a^olMjil + q)-^'^- (11) 

This is comparable to the duration of the plateau in the 
disruption rates shown in Figure [31 implying that binary 
evolution should not qualitatively change the plateau val- 
ues. A more sophisticated calculation that couples the re- 
sults of numerical scattering experiments with an evolv- 
ing binary will be the subject of a subsequent paper. 
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Fig. 3. — Numerical tidal disruption rates as a function of time for a hard MBH binary of mass ratio q, eccentricity e, and separation 
a = afi, embedded in an isothermal stellar cusp. The derivation assumes Mi = 10^ Mq, m» = 1 Mq, r* = Rq, and cr* = 100 kms^^. 
The short vertical lines mark the Kozai timescale for a» = a = aj^. 
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